Ultrafast Dynamic Metallization of Dielectric Nanofilms by Strong Single-Cycle 

Optical Fields 



Maxim Durach, 1 Anastasia Rusina, 1 Matthias F. Kling, 2 and Mark I. Stockman 3,2 

1 Department of Physics and Astronomy, Georgia State University, Atlanta, Georgia 30303, USA 
2 Max- Planck- Institut fiir Quantenoptik, Hans-Kopfermann-Strafle 1, D-85748 Garching, Germany 
3 Department of Physics and Astronomy, Georgia State University, Atlanta, Georgia 30303, USA 
E-mail: mstockman@gsu.edu; Homepage: \http:/ / www.phy- astr. gsu. edu/ 'stockman 

(Dated: January 19, 2013) 

We predict a dynamic metallization effect where an ultrafast (single-cycle) optical pulse with a 
< 1 V/A field causes plasmonic metal-like behavior of a dielectric film with a few-nm thickness. 
This manifests itself in plasmonic oscillations of polarization and a significant population of the 
conduction band evolving on a ~ 1 fs time scale. These phenomena are due a combination of both 
adiabatic (reversible) and diabatic (for practical purposes irreversible) pathways. 
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Latest advances in the ultrafast optics have recently at- 
tracted a great deal of attention. Ultrashort pulses have 
been successfully employed for monitoring and manip- 
ulation of electronic processes in atomic and molecular 
structures!^ Significant efforts have been directed toward 
exploring the potential of ultrashort (~ 100 as to ~ 1 — 10 
fs in duration) pulses in application to condensed mat- 
ter dynamicipHS in particular, to plasmonic metal and 
dielectric nanostructuresPEH 

We have recently predicted that dielectric nanofilms 
subjected to strong but sufficiently slow (adiabatic) elec- 
tric fields undergo a reversible change resembling a quan- 
tum phase transition to a state that exhibits metallic 
optical properties.^ We have called this phenomenon 
metallization. The minimum duration of such an adia- 
batic field depends exponentially on the thickness of the 
nanofilm and is in the range from > 10 fs to ~ 10 ns for 
a film thickness from a few nm to ~ 10 nmP^ 

Both from the fundamental point of view and for ap- 
plications to ultrafast nanoelectronics, the metallization 
by much faster optical fields is of great interest. In this 
Letter, we theoretically predict a new effect that we call 
dynamic metallization, where a single-cycle optical pulse 
incident on a ~ 2 nm dielectric nanofilm with a normal 
polarization and a field < f V/A causes a population of 
the conduction band and metal-like plasmonic polariza- 
tion oscillations on the optical-period time scale. This 
effect is caused by both adiabatic (reversible) and di- 
abatic (dissipative) excitation pathways involving band 
anticrossings and adiabatic evolution between them. 

A comprehensive solution of the ultrafast electron dy- 
namics in strong optical fields would require many-body 
quantum kinetics, rendering this problem extremely com- 
plicated. To simplify it, we rely on the fact that the 
characteristic inelastic electron-scattering time is on the 
order or greater than the surface plasmon decay time r n , 
which is r„ > 10 fs for metals - see, e.g., Fig. 1 (a) in 
Ref. 1131 Using ultrashort excitation pulses with duration 
t <C T n , we avoid any significant effect of the electron 
inelastic scattering. This allows us to treat the electron 



dynamics as Hamiltonian. The evolution of the system in 
this case is convenient to describe by the density matrix 

p(r',v ) t)= ]T*i(r',t)**(r,t) , (1) 

where if denotes the Fermi-surface state, i.e., the high- 
est occupied state for the zero- field Hamiltonian Hq, and 
^i(r',t) are the one-electron wave functions. These sat- 
isfy the Schrodinger equation ift?S>i = H(t)9i, where H(t) 
is the Hamiltonian depending on time t due to the optical 
field, and the dot denotes the derivative over t. 

Consider a thin nanofilm where the energy bands are 
split into subbands due to the quantum confinement in 
the direction normal to the film plane. We assume that, 
due to the material symmetry, the electron wave func- 
tion can be factorized into normal and parallel to the 
film. Assuming a normal optical electric field £ = £(t), 
the one-particle Hamiltonian of the transverse motion is 
H {£ ) = Hq — £d, where d is the dipole operator. 

Consider the adiabatic basis of states ipi(£) that di- 
agonalize the instantaneous Hamiltonian, H(£)ipi(£) = 
Ei(£)ipi(£), where Ei(£) are the adiabatic energies. We 
employ the Kronig-Penney model for an insulator with 
Eg = 4.8 eV band gap at the zero field (simulating di- 
amond) whose adiabatic subband energies Ei(£) of the 
valence (red) and conduction (blue) bands are shown in 
Fig. [11(a) as functions of the applied field £. All calcula- 
tions are done for a I = 2 nm thickness film. 

We expand wave functions in the adiabatic ba- 

sis, 9i(t) = Ej exp[-z^(i)]af (t)Vj (£(*)), where af {t) 
are the expansion coefficients with the initial condition 
0^(0) = Sij, and ^(t) = -(1/K) J Ej (£(t)) dt is the adi- 
abatic phase. Then the Schrodinger equation becomes 

af = - £ e jfe exp [-i<p Jk (t)] af , (2) 

Q jk = -£d jk (£) /E jk {£) , (3) 

where the adiabatic dipole matrix elements, tran- 
sition energies, and relative phases are dj k (£) — 
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FIG. 1: (a) Energy spectrum of the nanofilm as a function of 
the adiabatically applied electric field. The occupied valence 
subbands are shown in red, the empty conduction subbands 
are in blue, (b) The diabatic coupling matrix element Oik 
between band-edge subbands [see Eq. (J3J) ] for different pulse 
amplitudes So, as indicated on the panel. 
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(M£)Wj(£)) and E jk {£) = Ej(£) - E k (£), Vjk {t) = 
tpj (t) — ipk (t) , correspondingly. 

Under the adiabatic conditions^, a strong electric field 
causes the band gap E g to decrease. The valence and 
conduction bands experience anticrossing at the metal- 
lization field £ m = 0.75 V/A with the anticrossing gap 
AE m — 0.45 eV - see Fig. [T] (a). If the field is increased 
adiabatically above £ m , the electrons are adiabatically 
(reversibly) transferred to the conduction band (hole) 
states and in space across the film. This is the met- 
allization transition where the optical properties of the 
nanofilm resemble those of a plasmonic metalP^ If the 
field is adiabatically switched off, the system returns to 
its ground state. The condition of the adiabaticity is ev- 
ident from Eqs. |2])-([3| and is t v h/AE m , where t p 
is the time needea for the field £(t) to pass through the 
anticrossing, in full agreement with Ref. 

In the opposite case of a fast diabatic passage of the 
anticrossing, the Schrodinger equation ([2| for the valence 
and conduction band-edge subbands, v and c, can be in- 
tegrated yielding the population of the conduction band 



n c (t) = sin 2 8 uc [£(t)] n/n sb 



(4) 



where n is the electron density and n s b is the number 
of the occupied subbands, n s b = 9 in the present model. 
Such rapid fields, in contrast to the adiabatic easeful do 
not induce the metal-like polarization: there is no spatial 
population transfer across the nanofilm. 

We consider the electron dynamics of a nanofilm sub- 
jected to an ultrafast field where both the adiabatic and 
diabatic processes contribute. The fastest dynamics is 
driven by single-cycle light pulses with duration r ~ 1 fs, 
which have recently been achieved.^ Here, we model a 
single-cycle pulse £ (t) by a waveform^ 



£{t) =£ Q e- u2 (l-2u 2 ) 



t/r 



(5) 



where the amplitude is £q, and the pulse duration is r. 
The pulse integral is zero, £{t)dt = 0, as should be. 

For such a pulse, in Fig. Jl] (b) we display the diabatic 
coupling matrix element Q ik between the valence and 
conduction band-edge subbands. Note that G)^ oc 1/r. 



FIG. 2: Polarization V and conduction-band population n c as 
functions of time t for various So and r. Normalized pulse field 
S(t)/So is shown by blue line. Normalized population n c (t)/n 
(scaled x 100) is displayed by the red curve; the same in DA is 
shown by the dashed gray curve. Relative polarization V /So 
is displayed by the green line; the same in AA is given by the 
dashed yellow line. The pulse length is r = 0.85 fs for (a)-(c) 
and r = 3.4 fs for (d). 



The peaks of the diabatic coupling element are at the 
adiabatic metallization points (band-edge anticrossings) , 
and they grow with the excitation-pulse amplitude. 

Now let us turn to the dynamics of the system excited 
by an ultrashort pulse with r = 0.85 fs (the mean fre- 
quency hto = 2H/t — 1.55 eV). This is illustrated in 
Figs. [2] (a)-(c) where we show conduction band popu- 
lation n c and polarization V — Tr{dp}/V , where V is 
the nanofilm's volume, as functions of time t for different 
pulse-field amplitudes £q. For comparison, we also show 
the excitation waveform and results obtained in the adi- 
abatic (AA) and diabatic (DA) approximations. 

As Fig. [2] (a) shows, for field £ significantly below 
the adiabatic metallization threshold, £ m = 0.75 V/A, 
the AA polarization^ and DA population Q follow the 
pulse. The computed polarization V (green curve) is 
close to the adiabatic case except for a small delay and 
low-amplitude oscillations on the pulse trailing edge with 
frequency « E g /h. This is due to the short duration of 
the pulse, which leaves at the end a partial coherence 
between the valence and conduction bands. The popu- 
lation n c (red curve) is small and dramatically retarded 
with respect to both the pulse and the DA curve, which 
is characteristic of the perturbative excitation. 

At the threshold of the adiabatic metallization, E ss 
£ m , as illustrated in Fig.[2](b), the conduction band pop- 
ulation n c dramatically increases. The calculated depen- 
dence n c (t) (red line) agrees well with the DA, except for 
the residual population after the pulse. The significant 
deviation of n c (t) from the DA (gray dash line) starts 
at the moment of the second anticrossing on the trailing 
pulse tail (t/r m 0.1) where the diabatic coupling peaks 
- cf. Fig. [l] (d). This adiabaticity violation causes the 
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significant residual population (t 3> t) and is 

dependent on the adiabatic phase ip as will be discussed 
below in conjunction with Fig. |4j 

In Fig.[2](b), the polarization V(t) is retarded by an al- 
most quarter pulse length (sa 7r/2 in phase) with respect 
to the driving pulse, which implies a strong absorption. 
There also coherent oscillations after the end of the exci- 
tation p ulse. All this is characteristic of plasmonic metal 
systems I^HSJ The polarization oscillations exhibit beat- 
ings between the frequency of the interband and much 
slower intraband transitions. The latter are caused by 
the pulse imprinting its frequency by polarizing the hot 
carriers in the conduction band. We call this effect the 
dynamic metallization. It is an ultrafast and dissipative 
strong-field transition to a plasmonic metal-like behavior. 

A similar phase delay between the excitation field and 
the polarization oscillations has been computed and at- 
tributed to the appearance of free electrons in the time- 
dependent density-functional theory^ of breakdown in 
bulk dielectrics subjected to high optical fields. Note 
that such a breakdown for quasi-stationary fields was in- 
troduced by ZenerP^ Importantly, the present dynamic 
metallization in thin films is fundamentally different. It 
is based on the adiabatic contribution to polarization, de- 
pends critically on the film thickness, and occurs at much 
lower intensities: our field £ < 1 V/A corresponds to the 
radiation intensity / < 3 x 10 13 W/cm 2 , in contrast to 
I ~ 10 15 W/cm 2 in Ref.QH 

For the 0.85-fs pulse with amplitude £ = 0.94 V/A, 
which is significantly greater than the adiabatic 
metallization-threshold field £ m = 0.75 V/A [Fig.|](c)], 
the dynamic metallization phenomena become even more 
developed. The magnitudes of population n c and po- 
larization V increase. The field time dependence n c (t) 
shows a pronounced saturation between the metalliza- 
tion (anticrossing) points at t/r s» ±0.2. The residual 
population forms due to the adiabaticity violation at the 
anticrossing at t/r « 0.2 and is relatively large because 
of the large diabatic coupling at this instance - cf. the 
corresponding (blue) curve in Fig.[l](d). The polariza- 
tion shows a pronounced plasmonic metal-like behavior: 
an approximately quarter-oscillation delay with respect 
to the excitation pulse and the oscillations with a lower 
frequency, which is close to the pulse mean frequency u>o ■ 

The excitation dynamics for a longer pulse with r = 3.4 
fs is shown in Fig. [2] (d). The main difference from panel 
(c) is that the polarization V{t) peaks almost simultane- 
ously with the excitation pulse maximum, as characteris- 
tic of the adiabatic metallization - cf. Ref. |T2] Still there 
are the residual population n cr and oscillations of V{t) 
after the end of the excitation pulse, which imply non- 
adiabatic processes occurring at the level anticrossings. 

For the strongly-nonlinear and dispersive problem un- 
der consideration, a useful measure of the magnitude of 
the system's polarizability can reasonably be defined as 
the effective permittivity 
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FIG. 3: Effective permittivity £o as a function of pulse am- 
plitude So (a) and inverse pulse duration 1/t (b). The label 
AA denotes a result of the adiabatic approximation. 
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FIG. 4: Relative residual population of the conduction band 
n cr /n as a function of the pulse length r. The green curve is 
the result of the full numerical computation, the blue curve 
is the analytical result with two band-edge subbands, and 
the red with four. Arrows show the adiabatic phase ip. (a) 
£ = 0.94 V/A and (b) S = 1.2 V/A. 



where Vq is the maximum value of the polarization in 
the process. Note that the maximum of V{t) is gener- 
ally delayed in time with respect to that of £(t). For the 
shorter pulses in Figs. [2] (b)-(c), this delay is approxi- 
mately ~ 1/4 of the oscillation length (« ir/2 in phase), 
which implies a strong dissipation (a significant imagi- 
nary part of the complex permittivity). 

In Fig. [3] (a)-(b), we plot £0 as functions of the 
excitation-pulse amplitude £q and duration r. For very 
short pulses with r < 0.85 fs, e slowly increase with 
£0 due to contribution of perturbative nonlinear absorp- 
tion. The magnitude £0 ~ 20 is rather large because of 
the wide, high-frequency spectrum of the short pulses. 
Close to and above the adiabatic metallization thresh- 
old, £q > £ rn — 0.75 V/A, this effective permittivity 
dramatically increases for longer pulses with r > 2 fs, 
suggesting that it is dominated by the adiabatic metal- 
lization mechanism. In fact, for r = 5 fs, the dependence 
£o(£o) resembles that for the adiabatic permittivityI21 [cf. 
the blue and dashed blue lines in Fig. [3] (a)]. Note that 
the appreciable oscillations in the dependence of £0 on 
the pulse duration r seen for longer times in Fig. [3] (b) 
are due to the interference of the excitation amplitudes 
at the two anticrossings (at the leading and trailing edges 
of the pulse). These are analogous to the Ramsey fringes, 
as discussed below for the residual conduction-band pop- 
ulation n cr in conjunction with Fig. [4j 

In Fig. |4j we display the dependence of the residual 
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(after the pulse end) population n cr on the excitation 
pulse length r. A striking feature of this dependence is 
the presence of high-contrast oscillations. These have a 
very clear physical origin. In the adiabatic picture,^ 
when the pulse leading-edge field £{t) approaches the 
metallization threshold (causing the anticrossing of the 
adjacent valence- and conduction-band subbands), the 
valence electrons are shifted in space to one surface 
of the nanofilm in the direction of the field where the 
electrons occupy the quantum-bouncer states.^ When 
the field increases above this metallization threshold, 
the electron population is coherently transferred to the 
opposite surface.^ This creates polarization oscillating 
with the transition frequency between the valence- and 
and conduction-band edges, u) vc {t) = [E v (t) — E c (t)]/H, 
which then adiabatically evolves with time t. The phase 
accumulated by these oscillations between the time t m \ 
of the anticrossing passage at the leading pulse-edge and 
that at the trailing edge is ip — J^™ 2 ui vc (t)dt. 

If ip is such that the electrons at the moment t m i are 
shifted to the initial (in the direction of the maximum 
pulse field) surface of the nanofilm, then there is a large 
probability of their return back to the valence band, and 
the minimum of observed. Otherwise, the fringe 

maximum is reached. Thus these oscillations is analogous 
to the well-known Ramsey fringes. As indicated in Fig. 
[4j the adjacent minima and maxima of the n cr (r) fringes 
are indeed separated by the phase change At/? = n. As 
one can see from Fig. |4j these fringes are described ana- 
lytically reasonably well with two and very well with four 
band-edge subbbands taken into account. 

To conclude, in this Letter we have predicted a 



new effect: ultrafast dynamic metallization of dielec- 
tric nanofilms. A single-cycle ultrafast (duration ~ 1 
fs) optical pulse with the normal electric field of a < 
1 V/A amplitude incident on a dielectric nanofilm (here, 
a diamond-crystal film with thickness ~ 2 nm), induces 
a plasmonic metal-like dynamics that develops during an 
ultrashort period on the order of the pulse's duration. 
For pulses of 1 — 2 fs or longer, this dynamics is char- 
acterized by a large, metal-like polarization oscillating 
with optical frequencies. There is also a significant resid- 
ual population of the conduction band, which strongly 
depends upon and can be coherently controlled by the 
adiabatic phase <p accumulated between the two metal- 
ization instances at the leading and trailing edges of the 
excitation pulse. The polarization oscillations extend be- 
yond the pulse end and also depend on the accumulated 
adiabatic phase tp. Thus the dynamic metallization is 
due to the combination and mutual influence of both 
the rapid adiabatic (reversible) and diabatic (dissipative) 
mechanisms. This dynamic metallization effect can find 
applications in lightwave electronics,^ in particular, to 
create a field-effect transistor controlled by light's electric 
field with a ~ 100 THz bandwidth. 
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